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The fractional deviation of nose whistler group delay from the zero-temperature model has been 
calculated assuming a Maxwellian magnetosphere and a gyrofrequency electron density distribution; 
i.e., the electron density varies as the inverse cube with distance from the earth's center. The thermal 
correction to the zero-temperature group refractive index has been inserted into the nose-whistler 
jiroup delay integral to determine the modified group delay. Significant deviations from the zero- 
temperature group delay for ire(iu(Micies above the nose frequency liave been calculated. Since devi- 
ations from the zero-temperature dispersion of a few percent are readily discernible, nose wliistler 
data should provide at least an estimate of the upper bcmnd on magnetosphere tem{)erature. Twenty- 
three whistlers liave f)een analyzed and an upper bound on inagnetospliere leni|)<Mat ure of 2( 10') °K has 
been deteiniined. 



1. Introduction 

Whistlers are naturally occurring electromagnetic 
signals in the audio frequency range that result from 
the dispersion of lightning energy which has traveled 
through the outer ionosphere along the hues of force 
of the earth's magnetic field. The term nose whistler 
designates a particular class of whistlers that exhibits 
simultaneous rising and falling tones starting at a 
given frequency called the nose frequeticy. \\\v 
dispersion of the energy from the lightning source 
results from the action of free electrons in the presence 
of a circularly polarized electromagnetic wave (whose 
sense of rotation is the same as that of the free elec- 
trons in the ambient magnetic field). 



'Ttiis study was sponsored by the Stanford Research Institute. Some assistance was 
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A multipath whistler is a group of whistler traces 
originating from a single atmospheric that traverse 
different magnetosphere paths. A typical frequency- 
time characteristic of a multipath whistler is illustrated 
by the spectrogram of figure 1. The spectrogram 
shows the whistler amplitude, measured by the dis- 
play intensity, as a function of time (abscissa) and 
frequency (ordinate). The vertical lines in the spec- 
trogram indicate impulsive atmospherics. The 
horizontal lines indicate constant-frequency signals, 
usually power-line harmonics. The three strongest 
whistler traces of figure 1 are identified at the bottom 
of the spectrogram by the letters A, B, and C. 

Each of the three marked traces is in reality a group 
of several whistler-components, a fact most easily seen 
in Trace B. The nose frequency of the leading 
whistler of each of the three traces is indicated by the 
frequency indicated to the left of the spectrogram. 



EI 2 JULY 1963 1550=15 UT 




Figure 1. Typical nose whistler spectrogram. 
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The whistler traces are derived from a single atmos- 
pheric source some 1.78 sec prior to the arrival of the 
nose of Trace A. This source is not shown in the fig- 
ure, but its character is very much like that of many of 
the vertical lines seen in figure 1. 

The nose whistler upper cutoff frequency is rarely 
as much as 60 percent greater than the nose frequency. 
The upper cutoff frequency of nose whistlers is most 
often extremely sharp (Trace A and Trace C). The 
damping at cutoff is of the order of 20 dB for frequency 
changes about the cutoff frequency of less than 2 
percent. In many instances VLF emissions are trig- 
gered at the upper end of nose whistlers (Trace B); 
in a few cases, the upper cutoff is defined by a more 
gradual rate of change of attenuation of whistler 
energy with frequency. 

Recently, Scarf [1962] and Liemohn and Scarf 
[1962a, b] have examined the thermal attenuation of 
whistlers at the upper cutoff frequency in an effort to 
estimate the magnetosphere temperature. Assuming 
isotropic Maxwellian and Cauchy distributions of 
electron velocity, magnetosphere temperatures of 
105 o^ (Maxwellian) and 4 X 10^ °K (Cauchy) were 
estimated. In a later work, Liemohn and Scarf 
[1964] note that the mean thermal energy (tempera- 
ture) is probably lower than they had previously antici- 
pated. In the succeeding sections of this paper, the 
dispersion characteristics of nose whistlers will be 
used to set an upper bound on the magnetosphere 
temperature. 

There has been a considerable divergence in mag- 
netosphere temperature estimates. The earliest 
estimates, based upon the rate of escape of helium 
[Spitzer, 1948] placed the magnetosphere temperature 
between 1000 and 2000 °K. Similar temperatures for 
the lower magnetosphere have been proposed based 
upon satellite drag calculations [Harris and Jastrow, 
1959], incoherent backscatter experiments [VanZandt 
and Bowles, 1960], and solar Lyman-a line-width mea- 
surements [Bates and Patterson, 1961]. In contrast. 
Chapman [1957] estimates the temperature of the 
outer magnetosphere at 2 X 10^ °K. Chapman's 
estimate is predicated on an estended corona model 
derived from considering the heat outflow from 
the sun. Chapman further suggests a temperature 
gradient of 1 °K per km extending from the earth's 
surface to 5 earth radii. 



2. Theoretical Analysis 

The existence of nose whistlers has been explained 
by cold plasma theory; however, the zero-temperature 
analysis can be modified to include finite particle 
temperatures. The refractive index rip [Ginzburg, 
1961] of the "whistler" mode in a "temperate" plasma 
is given by 
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If the final bracket on the right-hand side of (1) is 
approximated by unity, the group refractive index, 
Ug^ becomes 

' df P''(fH-ff" 

a){alcff;pn{?>l2)fH + 2f] 
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and where 

f= wave frequency, 
/p = plasma frequency, 
/// = electron gyrofrequency, 

a= =root mean square electron speed, 

r= Electron temperature, 

/c^^Boltzmann's constant, 
m= electron mass, 

c= speed of light. 

The whistler time delay, t, can be determined at 
each frequency by an integration of the group refrac- 
tive index over the whistler path: 



r rigds 
is c ' 
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where ds is the elemental path length and the propa- 
gation is along the magnetic field. The whistler 
dispersion, D, is defined by: 



= tfl^ = Do^D, 
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Here Do defines the cold plasma dispersion, and 
Di the finite temperature perturbation to the zero- 
(1) temperature analysis. 
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The integrals may be evaluated subject to some 
assumptions on the magnetospheric whistler environ- 
nent and some changes of variable. Assume that: 

1) the earth is immersed in a dipole magnetic field; 

2) the magnetosphere has isothermal Maxwellian par- 
ade distribution; and (3) the magnetosphere electron 
density is proportional to the electron gyrofrequency, 
so that fp^fn^^^' The final assumption is con- 
sistent with the observations of Smith [1961] based 
upon whistler dispersion measurements below the 
nose frequency. 

Express the path length, S, in terms of the variable 
X, where x is the sine of the magnetic latitude, 6. 
If X is subsequently expressed in terms of the quantity 
y, where y is defined by 



1 



.y. 



(1-^y _fm 



(1 + 3:^2)1/2 / ' 



then the dispersion of a whistler propagating along a 
magnetic field line can be expressed by 
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Here ro is the radius of the earth, Oo is the latitude of 
the termination of the field line along the earth's 
surface, fnu is the gyrofrequency at the equator on the 
earth's surface, and m is the wave frequency,/, normal- 
ized to the minimum gyrofrequency along the path. 

The symbols y and h{y) are plotted in figure 2. 
For most latitudes of interest, yo is very close to unity 
(for ^0 > 55°, y > 0.99); consequently we shall assume 
yo= 1. The function h{y), too, is well behaved despite 
the singularity at y= 1. To evaluate the integrals of 
(5) and (6), we shall remove the factor h{y) from the 
integrals. The error introduced by these assumptions 
shall be considered below. 

The magnetosphere electron density is now assumed 
proportional to the electron gyrofrequency, so that 

where C is the constant of proportionality between 
density and gyrofrequency and is approximately 2(10^) 
when the electron density is 250 cm~'^ at 4 earth radii. 
The integrals of (5) and (6) can now be represented 
[T. F. Bell, private communication; Guthart, 1964] 



h { y ) 




Figure 2. Quantities y and h(y) as a function of magnetic lati- 
tude, 6. 



in terms of elliptic functions as shown below: 
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where E is the complete elliptic integral of the second 
kind and K is the elliptic quarter period [Milne- 
Thompson, 1950]. 

The normalized first-order disperson, (Z)i/Z)o), is 
expressible in terms of these elliptic functions and is 
given by 



Di^ 3{alcy'mC)F(m) 
Do 4/hi 



(10) 
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where 



F{m) = 



m(l - m) \l[{m) + {^ ml'[{m) 



E{m) 



Here F{m) is defined by the solutions of the elliptic 
integrals which have been evaluated using the well- 
known differential relations between the functions; 
F(m) is shown in figure 3 for 0.25 ^ ^n ^ 0.75. From 
the figure it is seen that F(m) is a rapidly varying func- 
tion of normalized frequency, m. 




0.8 



Figure 3. Elliptical integral weighting factor, F(m), as a function 
of normalized frequency. 



Equation (10) can now be used to measure the mag- 
netosphere temperature. Nose whistler data can be 
analyzed to determine whether there is a measurable 
increase in dispersion above the zero-temperature 
predictions. If there is a measurable increase in 
dispersion, then (10) measures the temperature 
directly. If no measurable increase is observed, (10) 
sets an upper bound on the temperature of the mag- 
netosphere, depending upon the accuracy to which 
the dispersion can be measured. 

The integrals of (7) through (9) have been evaluated, 
without approximations for yo and h{y), when m = 0.55 
and 771—0.65. The resulting solutions have been 
expressed in the form (DiIDq) and compared to the 
approximate solution of (10). The exact integral solu- 
tion is then seen to yield a value of F{m) that is 80 
percent of the value determined from the elliptic 
function approximation. 

A further examination of the integrands of the first- 
order dispersion integrals leads to the conclusion that 
at least 80 percent {m ^ 0.55) of the increased disper- 
sion introduced by temperature occurs within 10 deg 
(j= 0.352) of the geomagnetic equator. This justifies 
the statement that an isothermal magnetosphere is 
unnecessary and (10) measures the electron tempera- 
ture at the radial extremum of the magnetic field line. 



3. Dispersion Analysis of Nose Whistler 
Data 

When more than one nose whistler appears in a 
multipath whistler, the nose frequency of the succeed- 
ing nose whistlers almost always decreases in time, 
which is consistent with the observation that whistlers 
with lower nose frequencies traverse higher frequency 
field lines, and hence undergo more delay. The fact ^ 
that nose whistlers sometimes appear out of suc- 
cession in a whistler group has led to the deduction 
by Carpenter [1963] of a knee in the magnetospheric 
ionization density profile. 

The nose whistlers used in this analysis were ac- 
cumulated by Stanford University at Eights Station 
(EI), Antarctica, Byrd Station (BY), Antarctica, and 
Suffield Experimental Station (SES), Alberta. The 
techniques for recording and spectrographic data 
display are described in a Stanford Radioscience 
Laboratory report [Carpenter, I960]. A discussion of 
the detailed procedure for the evaluation of the spectro- 
graphic data leading to magnetosphere temperature 
estimates is found in Guthart [1964]. Briefly, the pro- 
cedure for data analysis is to scale the nose whistler 
spectrogram to determine: 

(1) The nose frequency, /n, 

(2) the upper cutoff frequency, /m, 

(3) the travel time to the nose from the whistler 
source, tn, 

(4) the travel time to the upper cutoff frequency from 
the whistler source, tm- 

The measured travel time at the upper cutoff fre- 
quency can be normalized to the travel time at the nose 
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frequency {tmlt,,)- The ratio can now be compared to 
the computed zero-temperature, dispersion curves 
of J. Angerami [private communication]. 

Angerami computed dispersion curves for an electron 
density model given by 

A^= electron density = Cxirojrf. 

In contrast, (10) is evaluated for a gyrofrequency elec- 
tron density model; i.e., 

I yV=C/,/ = C2(r«/r)^(l-h3 sin2 SY^K 

where tq is the radius of the earth, r is the radial dis- 
tance from the center of the earth, and 6 is the mag- 
netic latitude. Computations of normalized dispersion 
j using these two models for electron density yield 
i results differing by less than 1 percent {m < 0.7), 
I a neghgible error. More generally, the normalized 
I dispersion is not critically dependent on the assump- 
I tion that the magnetosphere electron density varies 
I as (ro/r)*^. If the magnetosphere electron density 
were, in fact, to vary as (to//)^ or (ro/r)'*, then, for the 
whistlers accumulated in this study (these whistlers 
shall be discussed shortly), the normalized dispersion 
would change by less than 2 percent from the (ro/r)*^ 
value. 

The comparison of the experimental and theoretical 
quantities results in a value for the increase in disper- 
sion introduced by the magnetospheric electron tem- 
perature; i.e., the comparison yields the quantity Di/Do, 
on the left-hand side of (10). The minimum value of 
DJDo is specified by the scaling accuracy, and is 
substituted into (10) when the experimental and the- 
oretical travel times are in agreement to the limit of 
the experimental accuracy. 

The scaled nose frequency and the travel time to the 
nose frequency are used in conjunction with the zero- 
temperature dispersion curves of Angerami to deter- 
mine the minimum gyrofrequency along the path and 
the constant of proportionality, C, relating the electron 
density and the electron gyrofrequency. The change 
in nose frequency introduced by a finite electron tem- 
perature is negligible for some typical nose whistler 
parameters [Guthart, 1964]. 

Equation (10) can now be solved for temperature. 
If a deviation in dispersion is observed, the inferred 
temperature is the electron temperature of the mag- 
netosphere. If no additional dispersion is observed, 
then we take the minimum observable dispersion from 
which we derive an upper bound on the magnetosphere 
electron temperature, since a higher temperature 
would have generated an observable dispersion. 

The accuracy of our estimate of magnetosphere 
temperature is critically dependent on the accuracy of 
our measurement of dispersion of the upper cutoff 
frequency. In reality, the measured quantity is the 
normalized dispersion, DmlD,,, i.e., the dispersion at 
the upper cutoff frequency normahzed to the dispersion 
at the nose frequency. Thus, DmIDn is given by 

/) t /f \i/2 

lJ m IjH ( JjH 



The greatest uncertainty in scahng dispersion is in 
the determination of the time of occurrence of the 
atmospheric from which the whistler originated, i.e., 
in determining tn. However, in the measurement of 
normalized dispersion, the desired quantity is tmltn, 
where 

tm *n~^L^ tm 

tn tn 

When /Mm = 0.65, A^m is only 15 percent of tn; 
consequently a 5-percent error in determining the 
travel time to the nose frequency, tn, generates an error 
of less than 1 percent in the scahng of (tmltn). A 5- 
percent determination of tn is readily achieved. 

The most common method of atmospheric (sferic) 
determination is based on the similarity of the traces 
of whistlers during a given 2-min period (at a given 
receiving site, the VLF recorder is usually turned on 
for 2 min of every hour), or sometimes over a period of 
many hours. Spectrograms of different whistlers are 
lined up so that the corresponding whistler traces are 
superimposed. The causative atmospherics corre- 
spond on the superimposed spectrograms. An alter- 
native method is to scale the normalized dispersion of 
the frequencies below the nose frequency and use the 
dispersion curves of Angerami to infer the sferic 
location. The application of this technique to several 
whistlers within a whistler group permits positive 
identification of the causative sferic. In this study, 
both the overlay and scaling techniques were applied 
so that the travel time to the nose could be determined 
to better than 3 percent. Carpenter [1960] describes 
several empirical methods of identifying whistler 
sources. 

The whistler nose frequency can be scaled to within 
3 percent and the upper cutoff frequency can be scaled 
to better than 1 percent. As a consequence of these 
scaling accuracies, it is estimated that the normalized 
dispersion, DmIDn, can be scaled to within 2 percent. 
The minimum observable value for Z)i/A), the quantity 
that specified the first-order change in dispersion 
introduced by finite electron temperature, is thus 
taken to be 0.02. 

Some 23 nose whistlers have been scaled and in 
no case has the measured dispersion deviated from 
that predicted by the cold plasma (electron tempera- 
ture =0) analysis; i.e., no temperature effects were 
observed. The whistlers used in this study are listed 
in table 1, including their respective properties. 
These whistlers have the following spectrum of 
parameters: 

(1) 0.57 ^rrim^ 0.75; 

(2) 2.65 kc/s ^///i ^ 13.3 kc/s; 6.9 ^ r/ro ^ 4; 

(3) 5.9(103)^C^ 35.8(10'^); 

where mm=fmlfH\ is the upper cutoff frequency nor- 
malized to the minimum gyrofrequency along the 
path, r/ro is the normalized distance from the earth's 
center, and C is the constant of proportionality relat- 
ing the electron density and ambient magnetic field 
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strength in the gyrofrequency model for electron 
density. The scaled whistlers include samples from 
the local summer and winter as well as over the 24-hr 
period. The data repeatability is such that within a 
2-min recording period, numerous whistlers generate 
similar temperature estimates, e.g., numbers 1 to 3, 
table 1. Indeed, over a period of hours much of the 
nose whistler data will be repeated, e.g., numbers 1 
and 6, table 1. Further, it should be noted that the 
data constituting this sample have been selected to 
demonstrate the existence of enhanced dispersion 
introduced by finite electron temperatures. Whis- 



tlers having an upper cutoff frequency below the nose 
frequency have not been included because, for a given 
magnetosphere temperature, these whistlers would 
give rise to less enhancement of the dispersion. 

The minimum observed upper bound on magneto- 
sphere electron temperature is 2(10^) °K (number 16, 
table 1). (For the whistlers of table 1, 
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Table 1. Experimental nose whistler data 








Date 


Time(UT) hrrmin/ 
sec into run 


Sta." 


///. , 


///I 


C 


p 


{tmltn) exp 
{tmltn)th 










kcis 




elm''^ cls'^ 


°K 




1 


2 July 1963 


15:50/10 


EI 


10.5 


0.64 


14.6(103) 


6.8(10^) 


0.995 


2 


2 July 1963 


15:50/10 


EI 


8.6 


.58 


12.1(103) 


1.3(10''^) 


1.003 


3 


2 July 1963 


15:50/10 


EI 


10.6 


.66 


14.6(103) 


5.8(10^) 


1.008 


4 


2 July 1963 


15:50/10 


EI 


8.6 


.65 


12.7(103) 


5.8(10^) 


0.991 


5 


2 July 1963 


15:50/102 


EI 


10.15 


.60 


12.7(103) 


1.2(105) 


1.018 


6 


2 July 1963 


14:50/54 


EI 


11.5 


.64 


14.8(103) 


6.9(10^) 


1.013 


7 


18 July 1963 


15:50/50 


BY 


5.9 


.59 


5.9(103) 


1.7(105) 


0.994 


8 


18 July 1963 


15:50/50 


BY 


2.94 


.60 


6.9(103) 


6.4(10") 


1.016 


9 


18 July 1963 


14:50/30 


BY 


2.65 


..57 


9.2(103) 


5.6(10") 


1.003 


10 


18 July 1963 


14:50/30 


BY 


5.25 


.57 


7.1(103) 


1.4(105) 


1.009 


11 


3 July 1963 


11:50/41 


EI 


8.05 


.65 


12.5(103) 


5.7(10") 


0.987 


12 


3 July 1963 


11:50/61 


EI 


8.85 


.58 


11.2(103) 


1.4(105) 


1.000 


13 


3 July 1963 


12:50/116 


EI 


7.6 


.59 


10.5(103) 


1.3(105) 


1.014 


14 


14 June 1962 


14:50/83 


ST 


11.8 


.67 


19.2(103) 


4.6(10") 


1.000 


15 


12 Feb. 1963 


17:35/29 


EI 


11.8 


.57 


14 (103) 


1.6(105) 


0.998 


16 


6 June 1963 


14:50/48 


EI 


9.85 


.75 


13.8(103) 


1.8(10") 


.974 


17 


6 June 1963 


14:50/48 


EI 


6.72 


.59 


16 (103) 


6.8(10") 


1.003 


18 


6 June 1963 


04:50/74 


EI 


6.75 


.57 


9.1(103) 


7.5(10") 


1.011 


19 


6 June 1963 


06:50/90 


EI 


5.1 


.59 


11.3(103) 


8.2(10") 


0.984 


20 


6 June 1963 


15:50/49 


EI 


13.3 


.67 


9.5(103) 


1.0(105) 


.983 


21 


6 June 1963 


16:50/86 


EI 


11.6 


.59 


13.5(103) 


1.4(105) 


.984 


22 


17 July 1962 


15:50/60 


BY 


1.6 


.59 


12.9(103) 


4.5(10") 


.987 


23 


11 July 1961 


08:50/37 


BY 


7.5 


.63 


7.5(103) 


1.0(105) 


.985 



"Stations El, Eights, Antarctica: BY, Byrd, Antarctica: ST. Stanford, Calif. 
'' T is an upper bound on temperature. 
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Figure 4. Error histogram for travel time to the upper cutoff fre- 
quency of nose whistlers. 



As a consequence, (2) is approximately valid.) Figure 
4 is a histogram depicting the observed error in the 
normalized travel time to the upper cutoff frequency. 
The number of nose whistlers exceeding the zero- 
temperature prediction was 10; the number whose 
time was less than that predicted by the zero tem- 
perature theory was 11. Only in one case was the 
difference in observed and predicted times greater 
than 2 percent. As a consequence, a statistical anal- 
ysis of a larger sample should improve the temperature 
estimate considerably. 

Some of the spectrograms studied in this analysis 
are shown in figure 5. The number identifying each 
whistler corresponds to the similarly numbered whis- 
tler of table 1. In analyzing the whistlers selected for 
temperature evaluation, care had to be exercised to 
exclude from the analysis VLF emissions triggered at 
the end of a whistler. (As an example, see number 15 
of fig. 5c.) The emission at the end of a whistler was 
nearly always characterized on the spectrogram by a 
change in slope (for example, see fig. 5e). However, 
in some cases, there was sufficient doubt as to the 
point of origin of the VLF emissions that the whistler 
trace had to be discarded (for example, see the upper 
branch of Trace B of fig. 1). 
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EI 3 JULY 1963 1150=41 UT 




(a) 



Figure 5a. Spectrograms of nose whistlers used in this analysis. 



EI 18 JULY 1962 1550=50 UT 




(b) 



Figure 5b. Spectrograms of nose whistlers used in this analysis. 
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EI 12 FEB 1963 1735= 29 UT 








^VLF EMISSION • 

; 









(c) 



Figure 5c. Spectrograms of nose whistlers used in this analysis. 



EI 6 JUNE 1963 0450=74 UT 




(d) 
Figure 5d. Spectrograms of nose whistlers used in this analysis. 
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EI 3 JULY 1963 1150^ 41 UT 




(e) 



Figure 5e. Spectrograms of nose whistlers used in this analysis. 



4. Concluding Remarks 

The increased dispersion introduced by finite elec- 
tron temperatures has been calculated. The results 
have been compared to experimentally observed whis- 
tlers, and the observed whistlers are seen to exhibit 
no temperature effects. As a consequence, an upper 
bound on the magnetosphere electron temperature 
has been set at 2(10^) °K, an order of magnitude below 
the previous bound set by Chapman. With the ac- 
cumulation of more data samples and the apphcation 
of statistical analysis, it should be possible to lower 
the above upper bound on temperature. 
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